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Abstract 

Motivated by multipath routing, we introduce a multi-connected variant of spanners. For that 
purpose we introduce the p-multipath cost between two nodes u and v as the minimum weight of a 
collection of p internally vertex-disjoint paths between u and v. Given a weighted graph G, a subgraph 
H is a p-multipath s-spanner if for all u, v, the p-multipath cost between u and v in H is at most s 
times the p-multipath cost in G. The s factor is called the stretch. 

Building upon recent results on fault-tolerant spanners, we show how to build p-multipath spanners 
of constant stretch and of 0(n 1+1/ ' fe ) edges 1 , for fixed parameters p and k, n being the number of nodes 
of the graph. Such spanners can be constructed by a distributed algorithm running in O(k) rounds. 

Additionally, we give an improved construction for the case p = k = 2. Our spanner H has 0(n 3 ^ 2 ) 
edges and the p-multipath cost in H between any two node is at most twice the corresponding one in 
G plus 0(W), W being the maximum edge weight. 
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1 Introduction 

It is well-known [ADD + 93] that, for each integer k ^ 1, every n-vertex weighted graph G has a subgraph 
H, called spanner, with 0(n 1+1 / k ) edges and such that for all pairs u, v of vertices of G, djj(u, v) ^ (2k — 
1) • dc{u, v). Here do(u, v) denotes the distance between u and v in G, i.e., the length of a minimum cost 
path joining u to v. In other words, there is a trade-off between the size of H and its stretch, defined here by 
the factor 2k— 1. Such trade-off has been extensively used in several contexts. For instance, this can be the 
first step for the design of a Distance Oracle, a compact data structure supporting approximate distance 
query while using sub-quadratic space [TZ05, BGSU08, BK06]. It is also a key ingredient for several 
distributed algorithms to quickly compute a sparse skeleton of a connected graph, namely a connected 
spanning subgraph with only 0(n) edges. This can be done by choosing k — O(logn). The target 
distributed algorithm can then be run on the remaining skeleton [BE10]. The skeleton construction can 
be done in 0(k) rounds, whereas computing a spanning tree requires diameter rounds in general. We refer 
the reader to [Pet07] for an overview on graph spanner constructions. 

However, it is also proved in [TZ05] that if G is directed, then it may have no sub-digraph H having 
o(n 2 ) edges and constant stretch, the stretch being defined analogously by the maximum ratio between the 
one-way distance from u to v in H and the one-way distance from u to v in G. Nevertheless, a size/stretch 
trade-off exists for the round-trip distance, defined as the sum of a minimum cost of a dipath from u to 
v, and a minimum cost dipath from v to u (see [CWOO, RTZ08]). Similar trade-offs exist if we consider 
the p-edge-disjoint multipath distance (in undirected graphs) for each p 1, that is the minimum sum of 
p edge-disjoint paths joining u and v, see [GGV10]. 
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1.1 Trade-offs for non-increasing graph metric 

More generally, we are interested in size/stretch trade-offs for graphs (or digraphs) for some non-increasing 
graph metric. A non-increasing graph metric S associates with each pair of vertices u, v some non-negative 
cost that can only decrease when adding edges. In other words, 8q{u,v) ^ 5h(u,v) for all vertices u, v 
and spanning subgraphs H of G. Moreover, if 5h(u,v) ^ a ■ 5g(u,v) + /3, then we say that H is an 
(a, /3)-spanner and that its stretch (w.r.t. the graph metric 5) is at most (a, (3). We simply say that H is 
an a-spanner if j3 = 0. The size of a spanner is the number of its edges. 

In the previous discussion we saw that every graph or digraph has a spanner H of size o(n 2 ) and with 
bounded stretch for graph metrics S such as round-trip, p-edge-disjoint multipath, and the usual graph 
distance. However, it does not hold for one-way distance. A fundamental task is to determine which graph 
metrics 8 support such size/stretch trade-off. We observe that the three former graph metrics cited above 
have the triangle inequality property, whereas the one-way metric does not. 

This paper deals with the construction of spanners for the vertex-disjoint multipath metric. A p- 
multipath between u and v is a subgraph composed of the union of p pairwise internally vertex-disjoint 
paths joining u and v. The cost of a p-multipath between u and v is the sum of the weight of the edges it 
contains. Given an undirected positively weighted graph G, define 5q(u,v) as the minimum cost of a p- 
multipath between u and v if it exists, and oo otherwise. A p-multipath s-spanner is a spanner H of G with 
stretch at most s w.r.t. the graph metric S p . In other words, for all vertices u, v of G, 5 p H (u 1 v) ^ s-6q(u, v), 
or 8 p H (u, v) < a ■ Sq(u, v) +(3 if s = (a, (3). It generalizes classical spanners as dc(u, v) — Sq(u, v) for p = 1. 

1.2 Motivations 

Our interest in the node-disjoint multipath graph metric stems from the need for multipath routing in 
networks. Using multiple paths between a pair of nodes is an obvious way to aggregate bandwidth. 
Additionally, a classical approach to quickly overcome link failures consists in pre-computing fail-over paths 
which are disjoint from primary paths [KKKM07, PSA05, NCD01]. Multipath routing can be used for 
traffic load balancing and for minimizing delays. It has been extensively studied in ad hoc networks for load 
balancing, fault-tolerance, higher aggregate bandwidth, diversity coding, minimizing energy consumption 
(see [MTG03] for a quick overview). Considering only a subset of links is a practical concern in link state 
routing in ad hoc networks [ JV09] . This raises the problem of computing spanners for the multipath graph 
metric, a first step towards constructing compact multipath routing schemes. 

1.3 Our contributions 

Our main contribution is to show that sparse p-multipath spanners of constant stretch do exist for each 
p ^ 1. Moreover, they can be constructed locally in a constant number of rounds. More precisely, we show 
that: 

1. Every weighted graph with n vertices has a p-multipath kp ■ 0(1 + p/k) 2k 1 -spanner of size 0(p 2 ■ 
rt 1+1//fc ), where k and p are integral parameters ^ 1. Moreover, such a multipath spanner can be 
constructed distributively in 0(k) rounds. 

2. For p = k = 2, we improve this construction whose stretch is 18. Our algorithm provides a 2- 
multipath (2, (9(W / ))-spanner of size 0(n 3 / 2 ) where W is the largest edge weight of the input graph. 

Distributed algorithms are given in the classical LOCAL model of computations (cf. [PelOO]), a.k.a. the 
free model [Lin92] . In this model nodes operate in synchronous discrete rounds (nodes are also assumed to 
wake up simultaneously). At each round, a node can send and/or receive messages of unbounded capacity 
to/from its neighbors and can perform any amount of local computations. Hence, each round costs one 
time unit. Also, nodes have unique identifiers that can be used for breaking symmetry. As long as we are 
concerned with running time (number of rounds) and not with the cost of communication, synchronous 
and asynchronous message passing models are equivalent. 

1.4 Overview 

Multipath spanners have some flavors of fault-tolerant spanners, notion introduced in [CLPR10] for general 
graphs. A subgraph H is an r-fault tolerant s-spanner of G if for any set F of at most r ^ faulty vertices, 
and for any pair u, v of vertices outside F, oIh\f{u, v) ^ s ■ oIq\f(u, v). 
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At first glance, r-fault tolerant spanners seem related to (r + l)-multipath spanners. (Note that both 
notions coincide to usual spanners if r = 0.) This is motivated by the fact that, if for an edge uv of G 
that is not in H, and if, for each set F of r vertices, u and v are connected in H \ F, then by Menger's 
Theorem H must contain some p-multipath between u and v. If the connectivity condition fulfills, there 
is no guarantee however on the cost of the p-multipath in H compared to the optimal one in G. Actually, 
as presented on Fig. 1, there are 1-fault tolerant s-spanners that are 2-multipath but with arbitrarily large 
stretch. 




Figure 1: A weighted graph G composed of a cycle of n + 1 vertices plus n — 1 extra edges, and a spanner 
H = G\ {uv}. Edge uv has weight 1, non-cycle edges have weight s, and cycle edges weight s/n so that 
cIh{u,v) — s. Removing any vertex z (fc {u,v} implies dQ\{ z y(u,v) = 1 and dn\{ z } ( w > v) — 2s(l — 1/rt). 
For other pairs of vertices x, y, dH\t z }(x,y)/dQ\r z \(x,y) < 2s. Thus, H is a 1-fault tolerant 2s-spanner. 
However $jj{u, v)/6q(u, v) > sn/s. Thus, H is a 2-multipath spanner with stretch at least n. 

Nevertheless, a relationship can be established between p-mutlipath spanners and some r-fault tolerant 
spanners. In fact, we prove in Section 2.4 that every r-fault tolerant s-spanner that is b-hop is a (r + 1)- 
multipath spanner with stretch bounded by a function of 6, r and s. Informally, a 6-hop spanner H must 
replace every edge uv of G not in H by a path simultaneously of low cost and composed of at most b edges. 
We observe that many classical spanner constructions (including the greedy one) do not provide bounded- 
hop spanners, although such spanners exist as proved in Section 2.1. Some variant presented in [CLPR10] 
of the Thorup-Zwick constructions [TZ05] are also bounded-hop (Section 2.2). Combining these specific 
spanners with the generic construction of fault tolerant spanners of [DK11], we show in Section 2.3 how 
to obtained a COCAC distributed algorithm for computing a p-mutlipath spanner of bounded stretch. 
A maybe surprising fact is that the number of rounds is independent of p and n. We stress that the 
distributed algorithm that we obtain has significantly better running time than the original one presented 
in [DK11] that was Q(p 3 \ogn). 

For instance for p = 2, our construction can produce a 2-multipath 18-spanner with (9(n 3 / 2 log 3 ' 2 n) 
edges. For this particular case we improve the general construction in Section 3 with a completely different 
approach providing a low multiplicative stretch, namely 2, at the cost of an additive term depending of 
the largest edge weight. 

We note that the graph metric S p does not respect the triangle inequality for p > 1. For p — 2, a cycle 
from utotu and a cycle from w to v does not imply the existence of a cycle from tt to u. The lack of 
this property introduces many complications for our second result. Basically, there are f2(n 2 ) pairs u, v of 
vertices, each one possibly defining a minimum cycle C u . v of cost 6q(u, v). If we want to create a spanner 
H with o(rt 2 ) edges, we cannot keep C u . v for all pairs u,v. Selecting some vertex w as pivot for going 
from u to v is usually a solution of save edges (in particular at least one between u and v) . One pivot can 
indeed serve for many other pairs. However, without the triangle inequality, C UtW and C w . v do not give 
any cost guarantee on 5 2 H {u, v). 

2 Main Construction 

In this section, we prove the following result: 

Theorem 1 Let G be a weighted graph with n vertices, and p,k be integral parameters ^ 1. Then, G 
has a p-multipath kp ■ 0(1 + p/k) " -spanner of size 0{kp 2 - l l k n l+l l k \og 2 - 1/k n) that can be constructed 
w.h.p. by a randomized distributed algorithm in O(k) rounds. 

Theorem 1 is proved by combining several constructions presented now. 



3 



2.1 Spanners with few hops 

An s-spanner H of a weighted graph G is b-hop if for every edge uv of G, there is a path in H between 
u and v composed of at most b edges and of cost at most s ■ iv(uv) (where ui(uv) denotes the cost of edge 
uv). An s-hop spanner is simply an s-hop s-spanner. 

If G is unweighted (or the edge-cost weights are uniform) , the concepts of s-hop spanner and s-spanner 
coincide. However, not all s-spanners are s-hop. In particular, the (2k — l)-spanners produced by the 
greedy 2 algorithm [ADD + 93] are not. 

For instance, consider a weighted cycle of n + 1 vertices and any stretch s such that 1 < s < n. All 
edges of the cycle have unit weight, but one, say the edge uv, which has weight u)(uv) — n/s. Note that 
da( u i v ) = ^(™) > 1. The greedy algorithm adds the n unit cost edges but the edge uv to H because 
dn{u,v) — n s • ui(uv) (recall that uv is added only if dn(u,v) > s ■ dc(u,v)). Therefore, H is an 
s-spanner but it is only an n-hop spanner. 

However, we have: 

Proposition 1 For each integer k ^ 1, every weighted graph with n vertices has a (2k — l)-hop spanner 
with less than n 1+1 ^ k edges. 

Proof. Consider a weighted graph G with edge-cost function u>. We construct the willing spanner H of 
G thanks to the following algorithm which can be seen as the dual of the classical greedy algorithm, till a 
variant of Kruskal's algorithm: 

(1) Initialize H with V(H) := V(G) and E(H) := 0; 

(2) Visit all the edges of G in non-decreasing order of their weights, and add the edge uv to H only if 
every path between u and v in H has more than 2k — 1 edges. 

Consider an edge uv of G. If uv is not in H then there must exist a path P in H from u to v such 
that P has at most 2k — 1 edges. We have d^(u, v) ^ w(P). Let e be an edge of P with maximum weight. 
We can bound oj(P) ^ (2k — 1) • u>(e). Since e has been considered before the edge uv, cu(e) uj(uv). It 
follows that uj(P) ^ (2k — 1) • uj(uv), and thus dn(u,v) (2k — 1) • oj(uv). Obviously, if uv belongs to H, 
dn(u,v) — uj(uv) ^ (2k — 1) • oj(uv) as well. Therefore, H is (2k — l)-hop. 

The fact that H is sparse comes from the fact that there is no cycle of length ^ 2k in H: whenever an 
edge is added to H, any path linking its endpoints has more than 2k — 1 edges, i.e., at least 2k. 

We observe that H is simple even if G is not. It has been proved in [AHL02] that every simple n- vertex 
m-edge graph where every cycle is of length at least 2k + 1 (i.e., of girth at least 2k + 1), must verify the 
Moore bound: 

fc-i 

n ^ l + dj^id- l) 4 > (d-l) k 

i=0 

where d = 2m/ n is the average degree of the graph. This implies that m < ^(n 1+1 ^ k + n) < n 1 + 1 / k . 

Therefore, H is a (2k — l)-hop spanner with at most n 1+1//fc edges. □ 

2.2 Distributed bounded hop spanners 

There are distributed constructions that provide s-hop spanners, at the cost of a small (poly-logarithmic 
in n) increase of the size of the spanner compared to Proposition 1. 

If we restrict our attention to deterministic algorithms, [DGPV08] provides for unweighted graphs a 
(2k — l)-hop spanner of size 0(kn 1+1 / k ). It runs in 3k — 2 rounds without any prior knowledge on the 
graph, and optimally in k rounds if n is available at each vertex. 

Proposition 2 There is a distributed randomized algorithm that, for every weighted graph G with n ver- 
tices, computes w.h.p. a (2k — T)-hop spanner of 0(kn 1+1 ^ k log 1 ^ 1 ^ n) edges in O(k) rounds. 

Proof. The algorithm is a distributed version of the spanner algorithm used in [CLPR10], which is based 
on the sampling technique of [TZ05]. We make the observation that this algorithm can run in 0(k) rounds. 
Let us briefly recall the construction of [CLPR10, p. 3415]. 

2 For each edge uv in non-decreasing order of their weights, add it to the spanner if d[j(u,v) > s ■ do(u,v). 
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To each vertex w of G is associated a tree rooted at w spanning the cluster of id, a particular subset 
of vertices denoted by C(w). The construction of C(w) is a refinement over the one given in [TZ05]. 
The main difference is that the clusters' depth is no more than k edges. The spanner is composed of the 
union of all such cluster spanning trees. The total number of edges is 0(/cn 1+1 / fe log 1_1/fc n). It is proved 
in [CLPR10] that for every edge uv of G. there is a cluster C(w) containing u and v. The path of the tree 
from w to one of the end-point has at most k — 1 edges and cost ^ (k — 1) • ui[uv), and the path from w 
to the other end-point has at most k edges and cost ^ k ■ u>(uv). This is therefore a (2k — l)-hop spanner. 

The random sampling of [TZ05] can be done without any round of communications, each vertex 
randomly select a level independently of the other vertices. Once the sampling is performed, the clusters 
and the trees can be constructed in O(k) rounds as their the depth is at most k. □ 



2.3 Fault tolerant spanners 

The algorithm of [DK11] for constructing fault tolerant spanners is randomized and generic. It takes 
as inputs a weighted graph G with n vertices, a parameter r ^ 0, and any algorithm A computing an 
s-spanner of m(v) edges for any z/-vertex subgraph of G. With high probability, it constructs for G an 
r-fault tolerant s-spanner of size 0(r 3 • m(2n/r) ■ logn). It works as follows: 
Set H := 0, and repeat independently 0(r 3 log n) times: 

(1) Compute a set S of vertices built by selecting each vertex with probability 1 — l/(r + 1); 

(2) H := HUA(G\S). 

Then, they show that for every fault set F C V(G) of size at most r, and every edge uv, there exists 
with high probability a set S as computed in Step (1) for which u,v (£ S and F C S. As a consequence, 
routine A(G\S) provides a path between u and v in G \ S (and thus also in G \ F) of cost $C s ■ ui(uv). If 
uv lies on a shortest path of G \ F, then this cost is s ■ (1q\f(u, v). From their construction, we have: 

Proposition 3 If A is a distributed algorithm constructing an s-hop spanner in t rounds, then algo- 
rithm [DKllj provides a randomized distributed algorithm that in t rounds constructs w.h.p. an s-hop 
r-fault tolerant spanner of size 0(r 3 ■ m(2n/r) ■ logn). 

Proof. The resulting spanner H is s-hop since either the edge uv of G is also in H, or a path between u 
and v approximating uj(uv) exists in some s-hop spanner given by algorithm A. This path has no more 
than s edges and cost s • u)(uv). 

Observe that the algorithm [DK11] consists of running in parallel q = 0(r 3 log n) times independent 
runs of algorithm A on different subgraphs of G, each one using t rounds. Round i of all these q runs can 
be done into a single round of communication, so that the total number of rounds is bounded by t, not by 

q- 

More precisely, each vertex first selects a g-bit vector, each bit set with probability 1 — l/(r + 1), its 
jth bit indicating whether it participates to the jth run of A. Then, q instances of algorithm A are run in 
parallel simultaneously by all the vertices, and whenever the algorithms perform their ith communication 
round, a single message concatenating the q messages is sent. Upon reception, a vertex expands the q 
messages and run the jth instance of algorithm A only if the jth bit of its vector is set. 

The number of rounds is no more than t. □ 



2.4 From fault tolerant to multipath spanner 

Theorem 2 Let H be a s-hop (p — l)-fault tolerant spanner of a weighted graph G. Then, H is also a 
p-multipath ip(s,p)-spanner of G where ip(s,p) = sp ■ 0(1 +p/s) s and ip(Z,p) = 9p. 

To prove Theorem 2, we need the following intermediate result, assuming that H and G satisfy the 
statement of Theorem 2. 

Lemma 1 Let uv be an edge of G of weight L)(uv) that is not in H . Then, H contains a p-multipath 
connecting u to v of cost at most ip(s,p) ■ uj(uv) where ip(s,p) = sp ■ 0(1 + pj ' sf and tp(3,p) = 9p. 
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Proof. From Menger's Theorem, the number of pairwise vertex-disjoint paths between two non-adjacent 
vertices x and y equals the minimum number of vertices whose removal disconnects x and y. 

By definition of H , H\F contains a path Pp of at most s edges between u and v for each set F of at 
most p — 1 vertices (excluding u and v). This is because u and v are always connected in G \ F, precisely 
by a single edge path of cost u>(uv). Consider Pjj the subgraph of H composed of the union of all such Pp 
paths (so from u to v in H \ F - see Fig. 2 for an example with p — 2 and s = 5). 

Vertices u and v are non-adjacent in Pjj. Thus by Menger's Theorem, Pjj has to contain a p-multipath 
between u and v. Ideally, we would like to show that this multipath has low cost. Unfortunately, Menger's 
Theorem cannot help us in this task. 

Let k s (u,v) be the minimum number of vertices in Pjj whose deletion destroys all paths of at most s 
edges between u and v, and let /i s (u, v) denote the maximum number of internally vertex-disjoint paths of 
at most s edges between u and v. Obviously, k s (u, v) ^ ^ s (u, and equality holds by Menger's Theorem 
if s = n — 1. Equality does not hold in general as presented in Fig. 2. However, equality holds if s is the 
minimum number of edges of a path between u and v, and for s = 2, 3, 4 (cf. [LNLP78]). 




Figure 2: A subgraph Pjj constructed by adding paths between u and v with at most s = 5 edges and 
with p = 2. Removing any vertex leaves a path of at most 5 edges, so k§(u,v) > 1. However, there 
aren't two vertex-disjoint paths from u to v of at most 5 edges, so k$(u,v) > fi§(u,v). Observe that 
/i 6 (w,u) = k 5 (u,v) = 2. 

Since not every path of at most s edges between u and v is destroyed after removing p — 1 vertices in 
Ph, we have that k s (u, v) ^ p. Let us bound the total number of edges in a p- multipath Q of minimum 
size between u and v in Pjj. Let r be the least number such that fx r {u, v) ^ p subject to n s (u, v) p. The 
total number of edges in Q is therefore no more than pr. 

By construction of Pjj, each edge of Pjj comes from a path in H \ F of cost lu(Pf) ^ s ■ d,Q\p(u, v) ^ 
s ■ uj(uv). In particular, each edge of Q has weight at most s ■ cu(uv). Therefore, the cost of Q is lu(Q) ^ 
prs ■ uj(uv). 

It has been proved in [PT93] that r can be upper bounded by a function r{s,p) < ( p j^2 2 ) + ( P ^2 3 ) = 
0(1 +p/s) s for integers s,p, and r(3,p) = 3 since as seen earlier k^(u,v) — /j, 3 (u,v). It follows that H 
contains a p-multipath Q between u and v of cost uj(Q) ^ sp ■ 0(1 + p/s) s ■ uj(uv) as claimed. □ 



Proof of Theorem 2. Let x, y be any two vertices of a graph G with edge-cost function u>. We want to 
show 5 p H {x,y) ^ Lp{s,p) ■ 5^{x,y). If 5 1 Q{x,y) — oo, then we are done. So, assume that 5q(x,j/) = lo(Pg) 
for some minimum cost p-multipath Pq between x and y in G. Note that lo{Pq) = ^2 uv ^e(p b ) 

We construct a subgraph Pjj between x and y in H by adding: (1) all the edges of Pq that are in H; 
and (2) for each edge uv of Pq that is not in H , the p-multipath Q uv connecting u and v in 7J as defined 
by Lemma 1. 

The cost of Pjj is therefore: 

w(-Pff) = ^2 UJ ( UV ) = ^2 u(uv) + ^2 u(Q uv ) 

uv&E(P H ) \uv£E(P G )nE(H) J \uv£E(P G )\E(H) 

By Lemma 1, w(Q uv ) ^ ip(s,p) ■ u(uv). It follows that: 

w(P H ) <p(s,p)- ^2 = • w(P G ) = ¥?(s,p) • 8 G {x,y) 

uv&E(P G ) 

as tp(s,p) 1 and by definition of P G . 

Clearly, all edges of Pp are in _ff . Let us show now that Pjj contains a p-multipath between x and y. 
We first assume x and ?/ are non-adjacent in Pjj. By Menger's Theorem applied between x and y in P#, 
if the removal of every set of at most p — 1 vertices in Pp does not disconnect x and y, then Pjj has to 
contain a p- multipath between x and y. 
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Let S be any set of less than p—1 faults in G. Since Pq is a p-multipath. Pq contains at least one path 
between x and y avoiding S. Let's call this path Q. For each edge uv of Q not in H , Q uv is a p-multipath, 
so it contains one path avoiding S. Note that Q uv may intersect Q wz for different edges uv and wz of Q. 
If it is the case then there is a path in Q uv U Q wz from u to z (avoiding v and w). assuming that u, v, w, z 
are encountered in this order when traversing Q. Overall there must be a path connecting x to y and 
avoiding S in the subgraph (Q D H) U UmieQ\-tf Quv- By Menger's Theorem, Ph contains a p-multipath 
between x and y. 

If a: and y are adjacent in Ph, then we can subdivide the edge xy into the edges x z and by adding 
a new vertex z. Denote by P' H this new subgraph. Clearly, if P' H contains a j>-multipath between x and 
y, then Pjj too: a path using vertex z in P' H necessarily uses the edges xz and zy. Now, P' H contains a 
p-multipath by Menger's Theorem applied on P' H between x and y that are non-adjacent. 

We have therefore constructed a p-multipath between x and y in of cost at most 
u(Ph) ^ v( s ;P) ' ^G^y)- ^ follows that 5 H (x,y) ^ tp(s,p) ■ S G (x,y) as claimed. □ 

Theorem 1 is proved by applying Theorem 2 to the construction of Proposition 3, which is based on 
the distributed construction of s-hop spanners given by Proposition 2. Observe that the number of edges 
of the spanner is bounded by 0(kp 3 ■ m(2n/p) ■ logn) = 0(kp 2 ~ 1 / k n 1+1 / k log 2-1 ^ n). 

3 Bi-path Spanners 

In this section we concentrate our attention on the case p = 2, i.e., 2-multipath spanners or bi-path spanners 
for short. Observe that for p = k = 2 the stretch is <p(3, 2) = 18 using our first construction (cf. Theorems 1 
and 2). We provide in this section the following improvement on the stretch and on the number of edges. 

Theorem 3 Every weighted graph with n vertices and maximum edge-weight W has a 2-multipath 
(2, 0(W))- spanner of size 0(n 3 / 2 ) that can be constructed in 0(n 4 ) time. 

While the construction shown earlier was essentially working on edges, the approach taken here is more 
global. Moreover, this construction essentially yields an additive stretch whereas the previous one is only 
multiplicative. Note that a 2-multipath between two nodes u and v corresponds to an elementary cycle. 
We will thus focus on cycles in this section. 

An algorithm is presented in Section 3.1. Its running time and the size of the spanner are analyzed in 
Section 3.2, and the stretch in Section 3.3. 

3.1 Construction 

Classical spanner algorithms combines the use of trees, balls, and clusters. These standard structures are 
not suitable to the graph metric S 2 since, for instance, two nodes belonging to a ball centered in a single 
vertex can be in two different bi-components 3 and therefore be at an infinite cost from each other. We will 
adapt theses standard notions to structures centered on edges rather than vertices. 

Consider a weighted graph G and with an edge uv that is not a cut-edge 4 . Let us denote by G[uv] the 
bi-component of G containing uv, and by S H (uv, w) the minimum cost of a cycle in subgraph H passing 
through the edge uv and vertex w, if it exists and oo otherwise. 

We define a 2-path spanning tree of root uv as a minimal subgraph T of G such that every vertex w 
of £r[m>] belongs to a cycle of T containing uv. Such definition is motivated by the following important 
property (see Property 1 in Section 3.3): for all vertices a, b in G[uv]\{u, v}, 5"q{cl, b) ^ ^(uu, a)+6^(uv, b). 
This can be seen as a triangle inequality like property. 

If 5t(uv,w) = Sq(uv,w) for every vertex w of G[iti>], T is called a shortest 2-path spanning tree. An 
important point, proved in Lemma 2 in Section 3.2, is that such T always exists and contains 0(v) edges, 
v being the number of vertices of G[itu]. 

In the following we denote by B 2 ; {uv, r) = {•«; : Sq{uv, w) S% ''} and Bq{u,t) = {w : dc{u,w) < r) the 
2-ball (resp. 1-ball) of G centered at edge uv (at vertex u) and of radius r. We denote by Na(u) the set 
of neighbors of u in G. We denote by BFS(m, r) any shortest path spanning tree of root u and of depth r 
(not counting the edge weights). Finally, we denote by SPST 2 <3(uz;) any shortest 2-path spanning tree of 
root uv in G[uw]. 

3 A short for 2- vertex-connected components. 

4 A cut-edge is an edge that does not belong to a cycle. 
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Figure 3: Process by which X is transformed into X'. 

The spanner H is constructed with Algorithm 1 from any weighted graph G having n vertices and 
maximum edge weight W . Essentially, the main loop of the algorithm selects an edge uv from the current 
graph lying at the center of a dense bi-component, adds the spanner H shortest 2-path spanning tree 
rooted at uv, and then destroys the neigborhood of uv. 

F := G, H := (0,0); 

while 3uv e E(G), \B G (uv,4W) n (N G (u) U N G (v))\ > do 
H := H U SPST 2 F (uv) U BFS G (u, 2) U BFS G (>, 2); 
|_ G := G \ (B 2 G (uv, AW) fl (N G (u) U N G (v))) 

H := HUG 

Algorithm 1: Construction of H . 

3.2 Size analysis 

The proof of the spanner's size is done in two steps, thanks to the two next lemmas. 

First, Lemma 2 shows that the while loop does not add too much edges: a shortest 2-path spanning 
tree with linear size always exists. It is built upon the algorithm of Suurballe-Tarjan [ST84] for finding 
shortest pairs of edge-disjoint paths in weighted digraphs. 

Lemma 2 For every weighted graph G and for every non cut-edge uv of G, there is a shortest 2-path 
spanning tree of root uv having O(v) edges where v is the number of vertices of G[uv\. It can be computed 
in time 0(n 2 ) where n is the number of vertices of G. 

Proof. In the following, we call X = G[uv}. SPST 2 x (uv) will therefore be equal to SPST 2 G (w). 
Let v = \V(X)\ and /z = \E(X)\. 

We first show that we can reduce our problem to finding a one-to-all pair of edge-disjoint paths in a 
directed graph. In other words, let V be a procedure which yields a 2-(edge-disjoint)-tree rooted in a single 
vertex w in a directed graph X' . We show we can derive V which yields SPST 2 js s :(mi') from V. 

First, remark that the problem of finding SPST 2 x(w , i;) is equivalent to finding the same structure but 
rooted in a single vertex w where the edge uv is replaced by uw, wv, and the weights of each edge uw and 
wv being equal to half of uj(uv). 

We construct X' as follows: each undirected edge is replaced by two edges going in opposite direction 
and of same weight. Then each vertex o is replaced by two vertices a\ and 02 where every incoming edge 
arrives at a\ and every leaving edge leaves from a^. An edge going from a\ to ai is finally added. Fig. 3 
shows what happens to edges of X. 

Note that v' = \V(X')\ = 2- {y + 1) and // = \E{X')\ = 2 • + 1) + v + 1. 

The procedure V' proceeds as follows: 

1. uv is replaced by uw, wv. 

2. X' is constructed. 

3. V is called on X' , with the root vertex being W2- 

4. Every edge of type X2 — > j/i present in the result of V causes the addition of the edge xy to the result 
of V. 
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Two edge-disjoint paths in X' are vertex-disjoint in X. Indeed, as they cannot both use an edge of the 
type x\ — > x 2 they cannot share X\ or x 2 (except at the extremities) because the only way to leave (resp. 
arrive) from x\ (resp. to x%) is to use the edge x\ — > x 2 . So if we have two edge-disjoint paths in X' going 
from W2 to some x±, the reduction back to X will yield two disjoint paths from w to x, and then from the 
edge uv to x. 

The procedure V was devised by Suurballe-Tarjan in [ST84]. While not directly constructing the 2- 
(directed-edge-disjoint)-tree rooted in a single vertex w it can be extracted from their algorithm. Roughly 
speaking, the construction results of two shortest-path spanning trees computed with Dijkstra's algorithm. 
The 2 paths (from the source to every vertex v) are reconstructed via a specific procedure. This latter can 
be analyzed so that the number of edges used in the 2-(directed-edge-disjoint)-tree is at most 2(z/ — 1): 
the structure Suurballe-Tarjan constructed is such that all vertices, but the source, have two parents. 

Therefore the number of edges yielded by procedure V is at most 4 • v, which is 0(v). □ 

Secondly, Lemma 3 shows that the graph G remaining after the while loop has only 0(n 3 / 2 ) edges. For 
that, G is transformed as an unweighted graph (edge weights are set to one) and we apply Lemma 3 with 
k = 2. The result we present is actually more general and interesting in its own right. Indeed, it gives an 
alternative proof of the well-known fact that graphs with no cycles of length ^ 2k have 0(n 1+1 / k ) edges 
since B G (uv, 2k) = in that case. 

Lemma 3 Let G be an unweighted graph with n vertices, and k J? 1 be an integer. If for every edge uv of 
G, \B G (uv, 2k) n N G (u)\ ^ n 1 ^, then G has at most 2 ■ n 1+1 ' k edges. 

Proof. Consider Algorithm 2 applied to graph G. When the procedure terminates, all the vertices and 
edges of the graph have been removed. In the following, we count the number of edges removed by each 
step of the while loop, which in turn allows us to bound the number of total edges of G. 



for i := k — 1 to do 

while 3u, \B G (u,i)\ ^ rW fe do 
[ L G:=G\B G (u,t) 

Algorithm 2: Remove 1-balls. 

Let Xi denote the set of vertices u whose 1-ball B G (u, i) is removed during iteration i of the for loop. 
Let m(u) be the number of edges erased when removing B G (u,i). Note that as X^Suex I^g( m j*)I = n 
(the procedure removed all the vertices), and that |JSQ| • n l l k ^ n because each 1-ball is larger than 
n l / k . ' 

At each step, we argue that 

m{u) ^ (n 1/fe + 1) • \B G (u,i)\ + \N G (u,i + 1)| 

where N G (u, i + 1) is the set of vertices at exactly i + 1 hops from u. 

To this effect, let's consider a shortest path tree T rooted in u and spanning B G (u, i). 

The number of edges in T is bounded by \B G (u,i + 1)|, which can be decomposed in \B G (u,i)\ + 
|JV G (u,i + l)|. 

We can also bound the total number of non-tree edges as follows: for any x £ B G (u,i), let's consider 
B G (xy, 2k) n N G (x), where y is the parent of x in T. We know that the number of vertices in this 2-ball 
is less than n x / k because it is a property of G. But \B G (xy,2k) H N G (x)\ is also at least the number of 
non-tree edges attached to x: for an edge xz T, the paths from z towards the root u and from x towards 
the root until they reach common vertex are of length at most the radius of B G (u, i), which is i ^ k, and 
so there is a cycle of length at most 2k using the edges xz and xy. So the total number of non-tree edges 
is bounded by n x l k ■ \B G (u,i)\. 

The termination of the while loop during iteration i + 1 implies 

\B G {u,i + l)\ < 

or equivalently: 

|JV G («,i + l)| < - \B G (u,i)\ 
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Therefore we have 

m(u) < n 1 / fc |i? G (u, i)\ + 

And so 

m(G) = J2 m{u)<n 1/k J2 J2 \B G (u,i)\ +J2\ X *\ ' n ^ 

ueUiXi i ueXi i 

and as Yli |Ai| ■ rc*/* ^ we have 

|#(G)| <2-n 1+1 / fe . 

□ 

Combining these two lemmas we have: 
Lemma 4 Algorithm 1 creates a spanner of size 0(n 3 ^ 2 ) in time 0(n 4 ). 

Proof. Each step of the while loop adds 0(n) edges from Lemma 2, and as it removes at least ypa vertices 
from the graph this can continue at most y/n times. In total the while loop adds 0(n 3 / 2 ) edges to H. 

After the while loop, the graph G is left with every B G (uv,4W) D (Ng{u) U Ng(v)) smaller than yfn. 
If we change all edges weights to 1. it is obvious that every B G (uv,4) n (Ng(u) U Ng(v)) is also smaller 
than s/n. Then as B G (uv,4) n N G (u) is always smaller than B G (uv,i) n (N G (u) U Ng(v)) we can apply 
Lemma 3 for k = 2, and therefore bound the number of edges added in the last step of Algorithm 1. 

The total number of edges of H is 0(n 3 / 2 ). 

The costly steps of the algorithm are the search of suitable edges uv and the cost of construction of 
SPST 2 . 

The search of suitable edges is bounded by the number of edges as an edge e which is not suitable can 
be discarded for the next search: removing edges from the graph cannot improve B G (e,4W). Then for 
starting from one extremity of each edge a breadth first search of depth 3 must be computed, keeping only 
the vertices whose path in the search encounters the other extremity. The cost of the search is bounded 
by the number of edges of G. So in the end the search costs at most 0(n 4 ). 

The cost of building a SPST 2 is bounded by the running time of [ST84], which at worst costs 0(n 2 ) 
(the reduction is essentially in 0(m + n)). Since the loop is executed at most yfn times, the total cost is 
0(n 7 / 2 ). 

So the total running time is 0(n 4 ). □ 



3.3 Stretch analysis 

The proof for the stretch is done as follows: we consider a, b two vertices such that Sp(a, b) = I is finite (if 
it is infinite there is nothing to prove) . We need to prove that the spanner construction is such that at the 
end, Sj I (a, b) ^ 21 + 0(W) . To this effect, we define Pp — Pp U Pp as a cycle composed of two disjoint 
paths (Pp and Pp) going from a to b such that its weight sums to Sp(a, b). 

Proving the stretch amounts to show that there exists a cycle Pjj = Pjj U Pjy joining a and b in the 
final H, with cost at most 2£ + 0(W) . Observe that if the cycle Pf has all its edges in H then one 
candidate for Pjj is Pp and we are done. If not, then there is at least one 2-ball whose deletion provokes 
actual deletion of edges from Pp (that is edges of Pf missing in the final H ) . 

In the following, let uv be the root edge of the first 2-ball whose removal deletes edges from Pp (that 
is they are not added in H neither during the while loop nor the last step of the algorithm). Let Gi be 
the graph G just before the removal of B G (uv, 4W) D (Ng(u) U Ng(v)) , and Gi+\ the one just after. 

The rest of the discussion is done in Gi otherwise noted. 

The proof is done as follows: we first show in Lemma 5 that any endpoint of a deleted edge (of Pp) 
belongs to an elementary cycle comprising the edge uv and of cost at most 6W. We then show in Lemma 6 
that we can construct cycles using a and/or b passing through the edge uv, effectively bounding 5ij(uv,a) 
and Sjj(uv,b) due to the addition of the shortest 2-path spanning tree rooted at uv. Finally we show in 
Lemma 7 that the union of a cycle passing through uv and a and another one passing through uv and b 
contains an elementary cycle joining a to 6, its cost being at most the sums of the costs of the two original 
cycles. 

Lemma 5 Let e = wt be an edge of (Gi \ G;+i) \ H . Then in Gi both w and t are connected to uv by a 
cycle of cost at most 6W. 
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Figure 4: Proof of Lemma 5: a cycle of 6 hops exists between w and uv in Gi 
u 




Figure 5: Proof of Lemma 7: the two cases for the simple paths. 



Proof. Since e isn't present neither in Gj+i nor H , then at least one of its endpoints is in the vicinity of 
u or v in Gi. W.l.o.g. we can consider t be a neighbour of u in Gi. w is then at most two hops from uv in 
Gi . The rest of the discussion is done in Gi . 

We can first eliminate the case where w is a direct neighbor of v, as there is an obvious cycle of 4 hops: 
w^t^u^v^w . 

Consider now the BFS tree rooted at u that is added to H . As w is at two hops at most from u there 
is a path u — > x — > w in this tree (x is defined as the intermediate vertex of this path and may not exist) . 
As e was removed, it means that x is distinct from t. Furthermore, t was removed because it belonged to 
a B 2 (uv,4W), so there is an elementary cycle of at most 4 hops passing through t and the edge uv. 

Now we distinguish two CctSGS cLS illustrated by Fig. 4. 

If x is distinct (which is especially true when it does not exist) from an intermediate vertex between v 
and t in the cycle, then we can extract an elementary cycle of at most 6 hops passing through uv and w : 
w — > x — > u — > v — >— > t — > w. 

If x is the same as an intermediate vertex between v and t, then the cycle vsw—tt—t-u—tv—tx—tw. 

□ 

We now show that we can use this lemma to exhibit cycles going from a to uv and from b to uv. 

From the vertices belonging to both Bq.(uv,6W) and Pp we choose the ones which are the closest 
from a and b (we know that at least two of them exist because one edge was removed from Pp during 
step i of the loop). There are at maximum four of them (oi, a2, &i, ^2), one for each sub-path Pp and each 
extremity {a, b}. Note that each extremity is connected to the root edge by an elementary cycle of cost at 
most 6W. Two cases are possible (the placement of the vertices can be seen on Fig. 5, although the paths 
on it are from the proof of lemma 7) : 

Case 1: There are only two extremities (then they belong to the same subpath) and their cycles which 
connect them to uv do not intersect the second subpath (w.l.o.g we can suppose it is a\ and 61). 

Case 2: There are more than two extremities: either some edges of the second path were removed or one 
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Figure 6: The two cases for Q\ and Q^. The big dots represent the paths' unused portions. 

of the cycles going from one of the extremities a± or b\ to uv intersects the second path. 

We show next that we can bound 8"jj{uv,a) and 6ff(uv,b) with the help of the cycles connecting the 
endpoints and the path Pq. This is done with the two next lemmas. 

Lemma 6 For any two vertices joined to the same edge uv by elementary cycles there is a simple path of 
cost at most the sum of the cycles' costs and passing through the edge uv. 

Proof. Let us call w and t the two vertices. We will show there is a simple path going from w to t passing 
through uv. Let us call Q 1 the elementary cycle joining w to uv and Q 2 the one joining t to uv. If by 
following Q 1 to reach from w one of the endpoints of uv we do not encounter Q 2 , then the path from w 
to t is composed of the part of Q , then the edge uv, then the part of Q 2 which reaches t without passing 
through uv. If it is not possible, then there are intersection points between Q 1 and Q 2 . Let i be the closest 
intersection point from w. The path we are looking for is therefore w — > i using Q 1 then i — > uv — > t using 
the part of Q 2 which uses the edge uv (the other part would take us directly to 02 whithout using uv). 
This path is simple because Q 2 is an elementary cycle and it cannot cross Q 1 before i because of the way 
i is chosen. The two cases are shown on figure 6. 

□ 



Lemma 7 Let a, b be two vertices such that an elementary cycle of cost S 2 (a, b) has common vertices with 
some B 2 (uv,6W). Then 5 2 (a,uv) and 5 2 lb,uv) are bounded by <5 2 (a, 6) + 12W 

Proof. The lemma is independent of the graph, but for clarity it will be proved using the graph Gi and 
Pf- 

Recall that we distinguished two cases depending on whether Bq.(uv,6W) intersects only one path of 
Pp (either P 1 p or P 2 p) or both. Fig. 5 illustrates the proof of the two cases. 

Lemma 6 allows us to solve the first case : since there are no intersection on the second path, the cycle 
a — Y o>i — ^ uv — y b{ — y b — y a is simple. So there is a cycle in Gi joining uv, a and b of cost at most 
12W + 6%.(a,b). So 8%.{uv,a) is bounded by UW + 6% {a,b) and so is 5q. (uv, b). 

In the second case there are three or four extremities: 01, 61, 0,2 and 62, with possibly 02 and &2 being 
the same vertex. We can apply Lemma 6 twice: once between a\ and 02 and another time between b\ and 
62 • These create a simple cycle from a to uv passing by a\ and 02 and another one from b to uv passing 
by bi and b 2 . We know the cycles are simple because the vertices were chosen to be the closest from a or 
b. Note that 02 and 62 can be the same. So 

5%. (a, uv) < u(a -!• 01) + 12W + w(a 2 -> a) < S 2 (a, b) + 12W 

and the same for b. 

□ 
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Property 1 Let uv be a non cut-edge of G and T be any 2-path spanning tree rooted at uv. Then, for all 
vertices a, b in G[uv] \ {u, v}, <5g(a, b) ^ S^(uv, a) + S^(uv, b) — uj(uv). 

Proof. There is in T a cycle joining a to uv of cost S^(uv,a), and another one joining b to uv of cost 
5^{uv,b). Consider the subgraph P containing only the edges from these two cycles. The cost of P is 
uj(P) ^ 5^(uv, a) + S^(uv, b) — tu(uv) as edge uv is counted twice. It remains to show that P contains an 
elementary cycle between a and b. Note that since a ^ {u, v}, a has in P two vertex-disjoint paths leaving 
a and excluding edge uv: one is going to u, and one to v. Similarly for vertex b. 

W.l.o.g. we can assume that a and b are not adjacent in P. Otherwise we can subdivide edge ab 
to obtain a new subgraph P' . Clearly, if P' contains an elementary cycle between a and b, then P too. 
Consider that one vertex z, outside a and b, is removed in P. From the remark above, in P \ {z}, there 
must exists a path leaving a and joining some vertex w a & {u, v} \ {z} and one path leaving b and joining 
some vertex wt € {u, v} \ {z}. If w a = Wb, then a and b are connected in P \ {z}. If w a ^ Wb, then 
edge uv belongs to P \ {z} since in this case z ^ {u, v}, and thus a path connected a to b in P \ {z}. By 
Menger's Theorem, P contains a 2-multipath between a and b. □ 

Lemma 8 H is a 2-multipath (2, 24W)-spanner. 

Proof. If there is in F a path of cost 5 2 (a,b) such that every edge of it is in H, then there is nothing to 
prove. If there is some removed edge, then we can identify the loop order i which removed the first edge, 
and we can associate the graph Gi just before the deletion performed in the second step of the loop (so 
Pf still completely exist in Gi). By virtue of Lemma 5 we can identify some root-edge uv and we know 
that there are some vertices of Pf linked to uv by an elementary cycle of length at most 6W. Lemma 7 
can then be applied, and so in Gi, c^. (», uv) and SQ.(b, uv) are both bounded by 5q. (a, b) + 12W . As the 
loop's first step is to build a shortest 2-path spanning tree rooted in uv we know that in H 

S 2 H (a,uv) < 5 2 Gi {a,uv) ^ S%.(a,b) + 12W 

and the same for b. Property 1 can then be used in the 2-path spanning tree, to bound 5 2 H {a, b): 

5 2 H (a, b) < S%(a, uv) + 5 2 H (b, uv) < 2 • 6%. (a,b) + 24W 

Finally, as in Gi Pf still exists completely, we have that Sq, (a, b) = 5p (a, b), so 

6 2 H (a,b) ^2-6 2 F {a,b) + 24W 

a 

4 Conclusion 

We have introduced a natural generalization of spanner, the vertex-disjoint path spanners. We proved 
that there exists for multipath spanners a size-stretch trade-off similar to classical spanners. We also have 
presented a 0(k) round distributed algorithm to construct p-multipath kp ■ 0(1 +p/k) 2k ^spanners of 
size 0(p 2 n 1+1 / k ), showing that the problem is local: it does not require communication between distant 
vertices. 

Our construction is based on fault tolerant spanner. An interesting question is to know if better 
construction (in term of stretch) exists as suggested by our alternative construction for p = 2. 

The most challenging question is to explicitly construct the p vertex-disjoint paths in the p-multipath 
spanner. This is probably as hard as constructing efficient routing algorithm from sparse spanner. We stress 
that there is a significant difference between proving the existence of short routes in a graph (or subgraph) , 
and constructing and explicitly describing such short routes. For instance it is known (see [GS11]) that 
sparse spanners may exist whereas routing in the spanner can be difficult (in term of space memory and 
stretch of the routes) . 
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